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ABSTRACT 

We present Herschel far-infrared (FIR) observations of two sub-mm bright quasars at high redshift: SDSS Jl 148-1-5251 (z = 6.42) 
and BR 1202-0725 (z = 4.69) obtained with the PACS instrument. Both objects are detected in the PACS photometric bands. The 
Herschel measurements provide additional data points that constrain the FIR spectral energy distributions (SEDs) of both sources, 
and they emphasise a broad range of dust temperatures in these objects. For /l,es, < 20yum, the two SEDs are very similar to the 
average SEDs of quasars at low redshift. In the FIR, however, both quasars show excess emission compared to low-z QSO templates, 
most likely from cold dust powered by vigorous star formation in the QSO host galaxies. For SDSS J 1148+5251 we detect another 
object at 160 ;um with a distance of ~ 10" from the QSO. Although no physical connection between the quasar and this object can be 
shown with the available data, it could potentially confuse low-resolution measurements, thus resulting in an overestimate of the FIR 
luminosity of the z = 6.42 quasar. 
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1. Introduction 

The detection of large quantities of cold dust in high redshift 
(z > 5) quasars (e.g. iBertoldi et all 120031: iBeelen et al.1 l2006t 
IWang et al]|2008l) implies a substantial enrichment of the inter- 
stellar medium already during the first billion years after the Big 
Bang. If the far-infrared (FIR) emission of these objects is pow- 
ered by star-formation, then the luminosity and the temperature 
of the dust implies star formation rates of up to a few thousand 
solar masses per year, possibly indicating the rapid formation of 
early galactic bulges. Recent observations at mid-infrared (MIR) 
wavelengths with Spitzer demonstrated that the dust also shows 
an energetically important hot (T ~ lOOOK) compon ent in many 
high-redshift QSOs (Hines et al l 120061: iJiang et all [2006), with 
only few exceptions (Jiang et an i2010l) . The similarity of the 
spectral energy distributions (SEDs) to those of lower redshift 
AGN lead to the conclusion that the general structures charac- 
terising local AGN are already in place at z ~ 6. However, for a 
majority of the sources the (sub)mm and MIR observations have 
only measured the tail of the dust emission spectrum. In order 
to further constrain the properties of the dust in these objects, it 
is essential to sample the SED as completely as possible. These 
measurements can then be used to derive critical parameters such 
as the total infrared luminosity and the range of dust tempera- 
tures in the objects of interest. In addition, the MIR to FIR lu- 
minosity ratio may indicate the relative importance of warm and 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA 



hot dust (which is predominantly heated by the AGN) compared 
to colder dust, which is preferentially heated by star formation. 
The Herschel key project (iPilbratt et al.ll2010l) "The Dusty 
Young Universe" (PI K. Meisenheimer) aims to measure the 
FIR SEDs for all quasars with z > 5 that were known at the 
time of submission of the proposal (early 2007) with th e PACS 
(iPoglitsch et all I20T0I) and SPIRE (Gri ffin et all I20T0I) instru- 
ments. Together with existing NIR and MIR photometry, mea- 
surements in five Herschel bands will yield complete infrared 
SED coverage of more than two decades in wavelength. Thus, 
the properties of the dust and the possible interplay between 
black-hole growth and galaxy bulge formation (traced indirectly 
through the FIR emission) can be explored early in the evolution 
of the universe. 



2. Observations and data reduction 

During the Science Demonstration Phase (SDP) two of our ob- 
jects with good ancillary data were observed with PACS. No 
SPIRE point-source photometry of the sources was obtained 
during the SDP. The objects were selected as SDP targets be- 
cause they were known to be bright at (sub)mm wavelengths. 
This allows us to study the performance of Herschel/PACS at 
low flux levels, and the experience gained here will be benefi- 
cial in the further course of the project when Herschel will be 
the only source of FIR data for many object. The observations 
also demonstrate HerscheFs ability to detect dust emission at 
the highest redshifts. 

S PSS J1148+525 1. The quasar SDSS Jl 148-1-5251 at z = 
6.42 (lFanetal.ll2003l hereafter J1148-H5251) was observed in 
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Fig. 1. Infrared and (sub)mm images of Jl 148+5251. All images are 60" wide and north is up with east to the left. Another nearby 
source northwest of the QSO (see text) is clearly visible at longer wavelengths. 



PACS scan-map mode using two cross scans in the green 
(lOO^ um, FWHM ~ 677) an d in the red (160A(m, FWHM -11") 
band (iPoglitsch et alJ2010l) . yielding an effective on-source time 
of ~1150sec in 3810sec of total AOR length. We employed 
standard processing procedures and masked the source location 
during deglitching and high-pass filtering. Both scan directions 
were processed separately and were then combined into a single 
map. The fluxes quoted in this paper were determined using an 
aperture of 5'.'0 radius (to avoid contaminations from neighbour- 
ing sources) and the sky was measured betwee n 30" and 35". 
App ropriate aperture corrections were included dPoglitsch et al.l 
1201 On . Due to the faintness of the source the fluxes measured 
via aperture photometry in the final maps were quite sensitive 
to the parameters chosen during the processing (i.e. the width of 
the high-pass filter and the size of the masked regions). In addi- 
tion to the scan maps, the source was also observed in chop-nod 
mode (~2100sec on-source time). Despite the higher noise lev- 
els in the chop-nod observations, the source is securely detected 
at 100 jum (but not at 160 yum). The scan maps with high-pass fil- 
ter widths of 15 and 20 in the green and red bands, respectively, 
yield fluxes that agree with the chop-nod observations. 

BR 1202-0725. The z = 4.69 quasar BR 1202-0725 
(lMcMahonetal.111994 hereafter 1202-0725) was observed in 
chop-nod mode in the blue (70;um, FWHM ~ 5"5) and red band 
of the PACS instrument for a total of 2550 sec (AOR length) with 
an on-source time of ~1980sec. Standard processing steps for 
chop-nod observations were used in combination with the latest 
calibration information. We performed aperture photometry us- 
ing apertures with a radius of 7'.'0 and KJ.'O in the blue and red 
filters, respectively. The sky was measured at distances between 
30" and 35" away from the source. 

3. Results 

3.1. J1148+5251 

The QSO Jl 148+5251 at z = 6.42 is one of the highest red- 
shift quasars know n to date. It has been detected previo usly 
at several sub-mm dBeelen et alJl2006l : iRobson et alJl2004l) and 



mm (iBertoldi et al.l l2003t iRiechers et al.l 12009 *) wavelengths. 
Under the assumption that the heating of the cold dust is dom- 
inated b y young stars, the high FIR luminosity (Lfir ~ 2 x 
10'^ Lq, ' Beelen et al.ll2006l) translates into a star- formation rate 
of ~ 3000 Moyr"'. These high star-formation rates are corrob- 
orated by measureme nts of the [ C II] emission line at 158 //m 
(iMaiolino et al.-20 05: Walte r et al.2009) . Mid-infrared observa- 
tions (iHines et al.l2 006; Jian g et al.l2006l) show large amounts of 
hot dust near the ~ 3 x 10^ M q black hole which accretes close 
to its Eddington Hmit (IWiUo tt et al. 2003; Barth et al. 2003). 

Figure[T] summarises the available multi-wavelength ob- 
servations in this field. A careful inspection of the data for 
Jl 148+5251 shows a companion ~ 10" to the northwest of 
the QSO in the FIR (this is best seen in the 160;um filter). 
Indeed, a similar extension is seen in the SHARC II observations 
(Beelen et al. 2006) and in new MAMBO measurements (based 
on .Bertoldi et al.,,2003.) of the source. The confusing source is 
clearly identified in many bands presented here, and a prelimi- 
nary analysis of deep HST/ACS imaging of this field does not re- 
veal an obvious count erpart. Interferome tric maps at ~1 mm ob- 
tained at IRAM/PdBI (IWalter et al.l2009h do not show a counter- 
part for the companion source. However, the observations were 
centred on the QSO, which places the location of the companion 
just outside of the primary beam at that wavelength, thus result- 
ing in highly reduced sensitivity. Given the currently available 
data on this source it is clear that further studies will be needed 
to reveal its true nature. We note though that depending on the 
shape of its SED this additional source will aff'ect measurements 




' see also the PACS Photometer release note available on the 
Herschel Science Centre website htt p : //herschel ■ esac ■ esa ■ int/| 



Fig. 2. Herschel/PACS chop-nod images of BR 1202-0725. The 
images are 120" wide and north is up with east to the left. 
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Fig. 3. Spectral flux distributions {top) and spectral energy distributions {bottom) for Jl 148+5251 {left) and for 1202-0725 {right). 
The arrow indicates the level of the 1.4 GHz (observed frame) radio emission (Tab.[Tl). In the top plot the solid line shows a single 
temperature grey body fit to the data at /lobs > 160/im. In the bottom panels we show as dashed lines the same greybody and also 
the mean SED for optically-luminous SDSS QSOs (Ric hards et alJf2006l) . scaled to the 1450 A (rest frame) data point. The mean 
SED was extrapolated for wavelengths longer than 100 //m (rest frame) assuming Fy ~ v^ (dotted line). The solid line is the sum 
of the scaled SED and the fitted greybody. Both objects clearly show a large excess of cold dust emission {A > 30//m, rest frame) 
compared to the mean QSO template. For J1148+5251 we here also include a Spitzer/IRS spectrum, which we retrieved from the 
archive and processed in a standard manner. 



of J1148H-5251 done with large (>10") beams. Aperture pho- 
tometry of J1148-I-5251 yields a flux density at 100 and 160 ;um 
of 4.1 + 1.0 mJy and 6.3 + 2.0 mJy, respectively (see Table[T]for 
a summary of all flux measurements, including literature data). 
In the red band we removed the confusing source via PSF fitting 
prior to performing the photometry. 



3.2. BR 1202-0725 

BR 1202-0725 is a well studied quasar at z = 4.69 wifli 
stron g detections at sub-mm and mm wavelengths (isaak et al.' 
[1994 : Omont et al. 1996b; Benford et al. 1999; lono et al. 2006) 
and in CO l ine emission ( Ohta et al.lll996t lOmont et al.lll996at 
ICarilU et ail [2002) . At ~4" distance to the norfliwest of flie 
quasar this source shows a secondary component in the dust con- 
tinuum and in CO. A Lyo- extension is observed approximately 
2''/3 northwest of the quasar (j Petitiean et al. 1996; Hu et al. 
11996; Fontan a et al.lll998l;IOhvam a et al. 2004). We detected BR 
1202-0725 at 70 /im and 160 /im with fluxes of 15.0 ± 2.0 mJy 
and 39. 8 + 3.7 mJy, respectiv ely (Fig.[2]l, which is consistent with 
the upper limits from ISO (Leech et al."200l'). In their Spitzer 
observations at 24yum, Hines et al. (2006) see slightly resolved 
emission and perform a two-component PSF fit to isolate the 
emission attributable to the QSO. Although our PACS lOfim 
data have slightly higher resolution than the MIPS 24 /im ob- 
servations, we do not detect a secondary component securely, 
probably due to the lower S/N in the 70jum maps compared to 
the 24 //m image. Although an apparent elongation may be visi- 
ble for the red band in Fig.|2] a two-component fit does not give 



a significantly better fit to the source profile. As we cannot re- 
liably separate the two sources in the PACS observations, the 
quoted flux values refer to the sum of both components. 



4. Discussion and conclusion 

In Fig. [3] we show the spectral flux and energy distributions of 
Jl 148+5251 and 1202-0725. The SEDs of both high-redshift 
QSOs appear very similar in shape, and the strong emission 
in all infrared bands shows dust at a wide range of temper- 
atures. This is particularly emphasised by the Herschel pho- 
tometry, which fills the gap between the previously available 
MIR and sub-mm photometry. The shape of the SED also im- 
plies that large amounts of dust are present in the host galax- 
ies already at high redshift and that the dust may be distributed 
on a wide range of scales: NIR emission from h ot dust close 
to the nucleus (e.g. iHines et al.l 120061: [jiang et al.l 12006, 2010|), 
warm MIR dust (T ~ few hundred K) on intermediate scales 
(or partly shielded), and colder dust in the FIR possibly dis- 
tributed throughout the host galaxy. Using the new Herschel 
photometry to determine total infrared luminosities for these 
two objects (integrating the SEDs between 1 and 200 //m, rest 
frame) yields Lhr = 8.1 x 10'^ Lq for Jl 148+5251 and a total 
Lfir = 3.7 X lO'^Lo for 1202-0725. 

The inspection of the spectral energy distributions in Fig. [3] 
shows that in the UV/optical and NIR/MIR (see below), where 
the emission is dominated by the active nucleus, the SEDs of 
the hi gh-z sources match th e template SEDs from lower redshift 
AGN dRichards et alJ2006l) reasonably well. However, while the 
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low-z infrared SEDs show only one broad peak at short MIR 
wavelengths, we observe a second peak with considerable ex- 
cess emission at FIR wavelengths in our targets. Greybody fits 
to the data at A > 160;um yield a temperature of ~ 50 - 60 K for 
the FIR dust emission, which is likely powered by the vigorous 
star formation (se e also lBertoldi et al.ll2003l; iBeelen et al.ll2006t 
IWang et aDl2008h . These large contributions to the FIR emis- 
sion, which can presumably be attributed to ongoing star forma- 
tion suggesting rapid bulge build-up in the host galaxies at high 
redshift, are missing in most lower redshift AGN. A combination 
of the mean QSO SED and the single temperature greybody rep- 
resenting the contributions from star formation is able to match 
the observed photometry. We note however that both objects 
presented here were selected as SDP targets because of their 
known strong FIR/sub(mm) emission. In principle, they could 
represent a small fraction of strongly star-forming objects, while 
the majority of the high-z QSOs may lack such a powerful FIR 
component. An analogue for this situation m ay be found in the 
(mostly local) PG quasars (iHaas et alJl2003h : when compared to 
the high-z objects, many of the PG sources show a similar mis- 
match in the FIR as seen for the SDSS QSO template, but a small 
number of PG quasars reveal a comparable FIR excess as seen 
for the high-redshift QSOs. Alternatively, stronger FIR emission 
from bulge build-up via star formation may be more common 
at high z than at low z, and we plan to explore these questions 
once data for a greater number of high-redshift objects become 
available. 

Both QSOs show more flux at NIR and MIR wavelengths 
than we would expect from the local templates. For 1202-0725 
this may be understood to be due to the inclusion of the com- 
panion (which is at the sa me redshift as the quasar: Omont et aP 
ll996aLICarilli et al.ll2002h and which contributes ro ughly half of 
the fl ux in many infrared a nd (sub)m m bands (e.g. lOmontet alj 
[T996a: Hines et al. 2006; lo no et alJ l2006). For J 1 148-1-5251 the 
situation is somewhat different, because the QSO is clearly the 
dominating source in flux at e.g. 8 and 24yum (observed, Fig.[T]). 
However, the exceptionally luminous black hole in this QSO ac- 
cretes close to its Eddington limit, which could possibly lead to 
a larger fraction of dust being heated to high temperatures, re- 
sulting in increased NIR emission. 

For both objects a secondary component may contribute to 
the measured flux densities at FIR wavelengths, potentially lead- 
ing to an overestimate of the FIR luminosity of the QSO itself. 
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